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HIGHLIGHTS 


•  A  physical  electric  circuit  model  was  used  to  model  two  similar  Li-ion  batteries  with  different  cathodes. 

•  The  use  of  this  model  allowed  calculation  of  exchange  current  density,  diffusion  coefficient,  and  double  layer  capacitance. 

•  This  approach  allowed  us  to  determine  the  corresponding  electrode  that  represented  these  behaviors. 

•  This  approach  was  non-destructive  to  full-cell  batteries. 
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A  physical  electric  circuit  model  (PECM)  was  used  to  identify  several  electrochemical  processes  occurring 
in  two  commercial  Li-ion  batteries  of  different  cathode  materials  (LixFeP04  and  LixCo02)  via  electro¬ 
chemical  impedance  spectroscopy  (EIS).  Through  defining  these  electrochemical  processes  in  these  two 
cells,  it  was  determined  that  the  charge  transfer  resistance  (or  exchange  current  density)  observed  via  EIS 
was  due  to  the  cathodic  exchange  current  densities  in  both  the  LixFeP04  and  LixCo02  full  cells.  In  dis¬ 
cussing  the  ionic  diffusion  of  the  examined  cells,  the  anode  of  one  cell  and  the  cathode  of  the  other  were 
primarily  responsible  for  the  observed  diffusion  of  the  full  cells.  Lastly,  the  measured  double  layer 
capacitance  was  determined  to  be  represented  in  EIS  scans  by  the  anodes  of  both  full  cells.  The  diffusion 
coefficient  was  calculated  using  Fields  1st  Law  estimation,  and  from  this  coefficient,  the  particle  size  was 
calculated  and  evaluated  against  scanning  electron  microscopy  (SEM). 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium-ion  batteries  (Li-ion  batteries  or  LIBs)  have  been 
commercially  available  for  more  than  20  years.  From  work  begun  by 
Goodenough  et  al.  and  Sony  Corporation  [1,2  ,  many  cathode  and 
anode  materials  have  been  investigated  for  use  in  their  construc¬ 
tion.  In  current  production,  two  cathode  chemistries,  namely  Lix_ 
FeP04  and  LixCo02,  have  grown  in  popularity  while  hard  carbon 
(HC)  and  graphite  have  been  used  primarily  as  anodes.  These  ma¬ 
terials  are  well  researched  with  their  materials  well  documented. 

To  date,  extensive  work  has  been  carried  out  in  determining  the 
physical  parameters  and  performance  of  battery  materials  in  a 
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variety  of  ways.  Zhou  et  al.  [3  have  used  electrochemical  imped¬ 
ance  spectroscopy  (EIS)  and  X-ray  diffraction  (XRD)  to  provide 
insight  into  nanoporous  cathode  mechanisms  including  phase 
transitions  of  cathode  materials  during  charge  and  discharge.  Osiak 
et  al.  [4]  reported  on  the  structural  and  electrochemical  charac¬ 
teristics  of  Sn02,  In203,  and  Sn  doped  In203  nanoparticles  used  as  a 
battery  anode  using  scanning  electron  microscopy  (SEM)  and 
energy-dispersive  X-ray  spectroscopy  (EDX)  techniques.  Levi  and 
Aurbach  5]  and  Della'Era  and  Pasquali  [6]  have  also  used  inter¬ 
mittent  titration  techniques  (PITT)  and  EIS  to  measure  the  diffusion 
coefficient  of  battery  materials,  and  have  shown  that  EIS-  and  PITT- 
measured  diffusion  coefficients  are  similar. 

Electrochemical  impedance  spectroscopy  (EIS)  is  a  mature  and 
useful  tool  in  characterizing  electrochemical  systems.  To  increase 
usefulness,  equivalent  circuit  modeling  has  been  employed  to 
provide  physical  or  phenomenological  representations  of  observed 
electrochemical  behavior.  Examples  of  the  usefulness  of  this 
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technique  can  be  found  many  places  in  literature:  Mauracher  and 
Karden  [7,8]  were  able  to  apply  this  modeling  to  lead-acid  batteries 
to  determine  important  battery  parameters  with  a  high  degree  of 
accuracy.  Stephen  Buller  [9]  also  applied  this  method  for  automo¬ 
tive  applications.  Moss  et  al.  were  able  to  use  this  method  to 
determine  various  features  of  energy  storage  devices  including 
electrochemical  capacitors  [10]  and  Li-polymer  batteries  [11-13]. 

In  this  paper,  a  physical  electric  circuit  model  (PECM)  discussed 
previously  14]  was  applied  to  EIS  data  gathered  on  two  commer¬ 
cial  Li-ion  batteries  of  similar  dimension  and  structure  (except  for 
cathode  material)  to  determine  if  it  is  possible  to  define  and 
quantify  the  electrochemical  processes  inside  a  Li-ion  battery  using 
only  EIS  and  other  in-situ  processes.  The  term  PECM  is  used  here  to 
define  an  electric  circuit  derived  from  a  physics-based  phenome¬ 
nological  model,  not  simply  an  electrical  analog  circuit  15].  This 
work  shows  the  significance  of  PECM  as  a  preferential  method  to 
analyzing  battery  performance  and  highlights  the  ability  of  PECM 
and  EIS  to  determine  the  factors  affecting  cell  performance. 

2.  Experimental 

Both  cells  used  in  this  experiment  are  commercially  available 
RCR123A  cells  from  PowerizerR  with  details  listed  in  Table  1.  The 
LixFeP04-cathode  cell  and  the  LixCo02-cathode  cell  were  assigned 
the  names  of  LFP  and  LCO  respectively.  As  was  done  previously  14], 
both  cells  were  subjected  to  potentiostatic  EIS  measurements  using 
a  Solartron™  1280B  Frequency  Response  Analyzer  at  uniform  states 
of  charge  (every  10%  SOC).  The  frequency  range  for  this  test  ranged 
from  20  kHz  to  10  mHz  with  0  V  DC  offset  and  an  AC  magnitude  of 
10  mV. 

3.  Modeling 

Select  resulting  spectra  for  both  cells  are  presented  in  Fig.  1. 
From  previous  work  [14]  the  EIS  data  was  fit  in  the  frequency 
domain  using  electric  circuit  elements  as  shown  in  Fig.  2.  While 
Ciucci  et  al.  have  developed  a  novel  EIS  parameter  evaluation 
methodology  [16,17  ,  each  element  was  fit  in  Scribner  and  Associ¬ 
ates'  Zplot®  software  using  a  Calc-modulus  data  weighing,  complex 
fitting,  iterative  method.  The  returned  accuracy  of  each  element  per 
EIS  scan  was  below  20%.  Also  it  should  be  noted  that  successive  EIS 


Table  1 

Measurement  of  battery  electrodes. 


Parameter 

LFP 

LCO 

Unit 

Capacity 

450 

650 

mAh 

Electrode  anode  length 

581.34 

499.14 

mm 

Electrode  anode  width 

25.12 

26.06 

mm 

Electrode  anode  thickness 

0.12 

0.17 

mm 

Electrode  cathode  length 

547.83 

482.83 

mm 

Electrode  cathode  width 

24.15 

24.79 

mm 

Electrode  cathode  thickness 

0.14 

0.14 

mm 

Separator  length 

982.35 

977.52 

mm 

Separator  width 

26.3 

26.309 

mm 

Separator  thickness 

0.02 

0.02 

mm 

Cu  current  collector  thickness 

0.01 

0.01 

mm 

Al  current  collector  thickness 

0.01 

0.01 

mm 

Anode  half  cell  weight 

3.6295 

4.2275 

g 

Cathode  half  cell  weight 

5.1594 

5.7876 

g 

Anode  weight 

2.324 

2.9904 

g 

Cathode  weight 

4.8022 

5.4661 

g 

Cu  volume 

0.146 

0.1383 

3 

cm 

Al  volume 

0.1323 

0.119 

3 

cm 

Cu  weight 

1.3055 

1.2371 

g 

Al  weight 

0.3572 

0.3214 

g 

Anode  electrode  area 

292.06 

260.15 

cm2 

Cathode  electrode  area 

264.6 

239.38 

cm2 

Fig.  1.  EIS  result  at  various  SOCs  for  (a)  LFP  and  (b)  LCO  cells. 


scans  of  the  same  cell  at  the  same  SOC  yielded  a  nearly  identical 
impedance  spectra  (hence  the  lack  of  error  bars  in  Fig.  1). 

By  assigning  certain  circuit  elements  to  electrochemical  pro¬ 
cesses,  we  may  evaluate  the  EIS  spectra  of  each  battery  and 
determine  (for  the  full  cell)  the  processes  that  affect  cell  perfor¬ 
mance.  These  elements  and  their  corresponding  processes  are 
defined  as  follows: 

L  describes  the  inductance  observed  in  the  Nyquist  plots  which 
correspond  to  the  inductance  of  the  experimental  leads  and  of  the 
spiral  wound  nature  of  the  cells.  Rs  describes  the  equivalent  ohmic 
series  resistance  of  separator,  solution  and  contact  resistances 
[18,19  in  the  cells.  Elements  R\  through  R3  and  C\  through  C3 
describe  the  formation  of  solid  electrolyte  interface  (SEI)  which 
form  on  the  surface  of  the  anode  [7,20,21  during  initial  cycling  and 
which  grow  through  successive  cycling.  Three  RC  circuits  were 
chosen  in  this  model  to  fit  the  SEI  layer  due  to  the  semi-elliptical 
nature  of  the  SEI  region  which  could  be  modeled  by  several  RC 
circuits  with  overlapping  spectra.  Rct  refers  to  the  charge  transfer 
resistance,  a  non-ohmic  reaction  resistance  that  relates  the  acti¬ 
vation  overpotential  and  exchange  current  density  as  described  by 
the  Butler-Volmer  equation  for  the  transfer  of  charge  carriers 
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Fig.  2.  Physical  equivalent  circuit  model  used  to  describe  EIS  data.  Note:  Rw  is  the  summation  of  rw  and  Cw  is  the  summation  of  cw.  For  this  circuit  (1 )  represents  the  line  and  spiral 
wound  inductance  of  the  cell,  (2)  represents  the  resistance  of  the  solution  and  separator,  (3)  represents  the  SEI  layer,  (4)  represents  charge  transfer  kinetics,  (5)  represents  ionic 
diffusion  into  cathode  and  anode  matrices,  and  (6)  represents  the  intercalation  capacitance  of  Li-ions  in  the  host  matrices. 


through  the  solid/electrolyte  boundary  and  Cdi  describes  the 
accumulation  of  charge  carriers  at  the  previously  stated  boundary 
condition  and  the  resulting  double-layer  capacitance  that  forms  as 
modeled  by  Guoy-Chapman  [20,22].  The  variables  7w,  i?w»  and  p 
refer  to  the  diffusion  time,  Warburg  resistance,  and  ideality  factor 
(respectively)  used  to  define  the  Warburg  equation  (Zw)  for  semi¬ 
infinite  diffusion  which  describes  the  semi-infinite  diffusion  of 
charge  carriers  through  media  14,15,22-24  .  Lastly  Qnt  defines  the 
intercalation  capacitance  of  the  host  materials. 

To  suitably  model  these  parameters,  equations  fitting  each 
parameter  as  a  function  of  SOC  were  developed.  Fig.  3  depicts  these 
equations  as  they  fit  parameter  as  a  function  of  SOC.  Based  upon 
these  equations,  and  through  the  application  of  Ohm's  law  and 
Inverse  Fourier  Transform,  we  can  determine  the  terminal  voltage 


of  the  battery  as  depicted  in  Fig.  4  which  is  governed  by  the 
following  equations: 


Vt(SOC)  =  U(  SOC) 


Vrs(SOC)  +  V'sei(SOC)  +  Vct||dl(SOC) 


+  ^cint(SOC)  +  ^Warburg  (SOC) 


where  U  is  the  open  circuit  voltage  of  the  cell  at  a  given  SOC, 


Vrs  (SOC)  =  I  x  1^5  (SOC) 


Fig.  3.  Fitting  equations  for  LFP  (solid  line)  and  LCO  (dashed  line)  equivalent  circuit  parameters.  Note  that  squares  represent  the  measured  LCO  parameters  and  circles  represent  the 
measured  LFP  parameter  where  graph  (a)  is  the  element  Rs,  (b)  Rc t,  (c)  Cdi,  (d)  Rw  and  (e)  Tw. 
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Fig.  4.  Time-domain  transform  of  the  frequency-domain  model  of  Fig.  2.  Note  that  U  is  the  open  circuit  potential  of  the  cell  at  rest. 


3 

\Zsm(SOC)  =  ^/xRn(SOC) 

"= 1  (3) 

x  (l  -exp(Rn(S0C)  xCn(S0C))) 

Vct||dl(SOC)  =  /xRct(SOC) 

X  (1_eXPUt(SOC)xCdl(SOC))) 

WSOC) , ,  X  Ks(SOC)  ,  (l  -  exp  (Bs(SOC,  ^(SOC,)  ) 

(5) 

oo 

vw,rburg(S0C)=J2ixRTii(S0C) 

n=1  (6) 

x  0  “  exp(%„(S0C)  x  cf  (soc) )  ) 

where  RTn  and  CTn  are  defined  by: 

k1(SOC)=Rw(SOC)  (7) 


4.  Discussion 

Under  the  paradigm  that  the  PECM  can  approximate  the  pa¬ 
rameters  governing  certain  electrochemical  processes  in  each  cell 
well  (with  respect  to  the  full  cell),  this  discussion  will  attempt  to 
ascribe  equivalent  circuit  model  parameters  to  phenomenological 
processes  occurring  in  the  cell  and  will  show  that  these  recorded 
parameter  values  are  reasonable.  The  scope  of  this  work  will  cover 
the  charge  transfer  kinetics  occurring  in  the  Li-ion  batteries  (i.e. 
charge  transfer  resistance  and  electric  double  layer  capacitance) 
and  the  ionic  diffusion  of  Li-ions  in  the  electrode(s).  This  discussion 
does  not  cover  the  intercalation  capacitance,  solid  electrolyte 
interface  (SEI)  and  the  ohmic  solution  resistance  of  the  cells  that 
was  calculated  in  the  Modeling  section  of  this  paper.  In  the  case  of 
the  former  two,  there  was  not  an  established  methodology  to 
validate  our  parameters  against  that  is  within  our  ability  to  test.  In 
the  latter  case,  the  validation  of  the  solution  resistance  was  trivial  in 
this  case  since  the  electrolyte,  separator,  and  contacts  were  the 
same. 

For  the  remaining  circuit  parameters  the  following  observations 
could  be  made:  The  charge  transfer  resistance  (Fct)  as  defined  as  the 
speed  of  the  surface  reaction  that  facilitates  the  transference  of  ions 
across  an  electrode/electrolyte  boundary  [26  may  also  be  repre¬ 
sented  by  the  exchange  current  density  via: 


/<2(SOC) 


Kw(SOC) 

Tw(SOC) 


CtJ  SOC) 


k-i  (SOC) 
2/<2  (SOC)2 


RtA  SOC) 


8/<|(SOC) 

(2 n  -  l)2it2 


(10) 


where  RT  and  CT  are  elements  of  the  transmission  line  model 

1  n  1n 

representation  of  the  Warburg  diffusion  element  as  discussed  by 
Levi  and  Aurbach  [20,25]  and  Maraucher  and  Karden  [7].  The 
equivalent  transmission  line  model  for  the  Warburg  diffusion 
element  can  be  seen  in  the  insert  in  Fig.  4. 

Based  upon  this  transformation  the  PECMs  for  the  LCO  and  LFP 
batteries  were  simulated  in  MATLAB™  against  experimental  dis¬ 
charges  carried  out  on  each  cell  using  an  Arbin  BT2000.  Fig.  5a  and 
b  shows  the  experimental  and  simulated  discharge  curves  for  each 
cell  at  various  C-rates,  while  Fig.  6  depicts  the  effect  of  select 
voltage  drops  (Eqs.  (2)— (10))  during  a  1  C  discharge  simulation 
provide  better  understanding  of  the  parameter  effect  as  a  function 
of  terminal  voltage. 


*o 


RT 


nFR 


ct 


(11) 


where  R  represents  the  universal  gas  constant,  T  represents  the 
temperature  in  Kelvin,  n  is  the  valence,  F  is  Faraday's  constant,  and 
Rct'  is  the  charge  transfer  resistance  per  unit  area. 

Table  2  lists  the  charge  transfer  resistance  and  corresponding 
exchange  current  density  of  the  LCO  and  LFP  batteries  as  a  func¬ 
tion  of  SOC.  The  respective  averages  for  the  LCO  and  LFP  battery 
exchange  current  densities  are  8.75  x  1CT3  and  2.7  x  1CT3  A  cm-2, 
respectively.  From  the  data  recorded,  two  conclusions  may  be 
made:  First,  the  exchange  current  density  of  the  full-cell  LCO  and 
LFP  batteries  is  on  the  same  order  of  magnitude  as  reported 
[25,27]  half-cell  (LixCo02,  LixFeP04,  Lii_xCe)  exchange  current 
densities,  leading  to  the  conclusion  that  the  model  interpretation 
of  the  charge  transfer  resistance  is  accurate.  Second,  the  general 
trend  of  marginally  decreasing  exchange  current  density  during 
cell  charging  leads  one  to  conclude  that  the  exchange  current 
densities  are  at  least  partially  dominated  by  their  cathode 
electrodes. 

The  first  conclusion  may  be  explained  in  greater  detail:  Values 
reported  for  LixCo02  by  Garcia  et  al.  [26  ranged  from  2.5  x  10-3  to 
3  x  10-4  A  cm-2  underperformed  compared  to  calculated  results 
taken  from  EIS.  This  may  be  due  to  the  difference  between 
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Fig.  5.  Experimental  and  simulation  results  for  (a)  LFP  cell  at  various  C-rates  and  (b) 
LCO  cell  at  various  C-rates. 


experimental  and  commercially  manufactured  cell  construction, 
though  the  inherent  error  used  in  modeling  Rct  via  EIS  may  have 
also  contributed  to  this  discrepancy.  In  the  case  of  the  LixFeP04,  the 
observed  exchange  current  density  of  2.7  x  1CT3  A  cnrT2  is  higher 
than  reported  by  others  [31,32]  (8  x  10-5  A  cm-2)  which  may  be 
explained  again  by  the  difference  in  producing  cells  in  a  laboratory 
vs.  a  commercial  process. 

The  second  conclusion  may  be  explained  further  by  relating  the 
decreasing  concentration  of  Li-ions  at  the  cathode  to  a  decrease  in 
the  magnitude  of  cathodic  exchange  current  densities.  This  trend  is 
well  documented  [26-31  when  looking  at  half-cell  performances. 
From  this  we  may  infer  that  the  trends  we  observed  are  caused  by 
the  cathodic  exchange  current. 

To  model  the  electrochemical  process  of  diffusion,  the  decision 
was  made  to  use  the  Finite  Length  Warburg  (FLW)  element  which 
models  1-D  particle  diffusion  with  a  transmissive  boundary  con¬ 
dition.  While  the  diffusion  of  Li-ions  in  LIBs  is  a  vastly  more  com¬ 
plex  system,  especially  in  dealing  with  two-phase  systems  like 
LixFeP04  and  LixCo02,  previous  work  5,6,34]  has  provided  good 
agreement  between  EIS  techniques  and  other  methods  (chiefly 
PITT)  in  determining  diffusion  aspects  including  the  diffusion  co¬ 
efficient  in  LIBs. 

In  incorporating  the  FLW  into  our  electric  model,  the  Warburg 
coefficient  (Aw)  was  determined  via  direct  calculation  (Eq.  (12))  and 


Fig.  6.  Graphical  representation  of  Eqs.  (2)-(9)  for  (a)  LCO  full-cell  and  (b)  LFP  full  cell 
for  1  C  discharge  as  a  function  of  terminal  voltage. 

validated  via  the  graphical  method  by  which  the  EIS  diffusion  tail  is 
plotted  Re(Z)  v.  1  /(w)~1/2  or  Im(Z)  v.  1  /(w)~1/2  where  the  slope  of  the 
resulting  lines  is  the  Warburg  coefficient.  The  difference  between 
graphical  and  calculated  Warburg  coefficients  are  7.45%  and  11.46% 


Table  2 

Charge  transfer  resistance  and  exchange  current  density. 


soc 

LFP 

Rct (£2) 

LCO 

Rct  (O) 

LFP 

i0  (A  cm'2) 

LCO 
*0  (A 

cm  2) 

0.1 

0.029196 

0.010925 

3.06 

X 

103 

9.01 

X 

10'3 

0.2 

0.032104 

0.00828 

2.79 

X 

103 

1.19 

X 

10'2 

0.3 

0.032888 

0.011522 

2.72 

X 

10'3 

8.54 

X 

10  3 

0.4 

0.030225 

0.010368 

2.96 

X 

10'3 

9.49 

X 

10'3 

0.5 

0.032384 

0.012162 

2.76 

X 

10'3 

8.09 

X 

10'3 

0.6 

0.030461 

0.01163 

2.94 

X 

10'3 

8.46 

X 

10  3 

0.7 

0.030568 

0.011847 

2.93 

X 

10'3 

8.31 

X 

10  3 

0.8 

0.031983 

0.010134 

2.80 

X 

10'3 

9.71 

X 

10'3 

0.9 

0.03282 

0.011531 

2.73 

X 

10'3 

8.54 

X 

10'3 
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for  the  LCO  and  LFP,  full  cells  respectively.  This  difference  further 
shows  good  agreement  between  the  PECM  and  EIS  data. 


RT  _  i?vv 
n2XF2SVW  ~ 


(12) 


where  X  is  the  molar  concentration  of  the  species,  S  is  the  surface 
area  of  the  electrode,  D  is  the  diffusion  coefficient  of  the  species, 
and  and  7w  are  the  Warburg  resistance  and  Warburg  time 
constant  respectively. 

With  the  Warburg  coefficient  for  the  LCO  and  LFP  full  cells 
accurately  determined,  it  may  be  possible  to  determine  if  one 
electrode  dominates  the  diffusion  coefficient  of  a  full  cell.  Ye  et  al. 
[35]  has  shown  that  the  Warburg  resistance  element  (i?w)  has  an 
inversely  proportional  relationship  to  the  diffusion  coefficient  (DL 0 
in  half  cells.  Also,  numerous  sources  [36-41  report  a  relationship 
in  which  the  diffusion  coefficient  in  electrode  materials  are 
inversely  proportional  to  the  amount  of  lithium  intercalated  into 
the  electrodes  (i.e.  Du  decreases  as  x  in  LixFeP04,  LixC6,  or  LixCo02 
increases).  In  the  case  of  the  LCO  cell  (as  seen  in  Table  3)  the 
Warburg  resistance  increases  with  increasing  SOC.  From  the  re¬ 
lationships  discussed  above,  if  the  Warburg  resistance  is  increasing 
as  Li-ions  are  being  shuttled  into  the  anode  (ion  concentration  is 
increasing  and  diffusion  coefficient  is  decreasing)  then  the  EIS  data 
presented  is  indicative  of  anode  behavior.  Similarly  for  the  LFP  full¬ 
cell,  when  seeing  an  increase  in  Warburg  resistance  with 
decreasing  SOC  it  may  be  inferred  that  the  EIS  data  is  reflecting  the 
cathode.  Table  3  lists  the  calculated  diffusion  coefficient  as  a 
function  of  SOC,  i?w  as  a  function  of  SOC,  and  the  calculated  and 
graphically  determined  Warburg  coefficients  (Aw)  as  a  function  of 
SOC  for  the  LCO  and  LFP  cells. 


(,3) 

Based  upon  this  inference,  the  diffusion  coefficient  for  each  cell 
may  be  calculated  as  shown  in  Eq.  (13).  Based  upon  assumed  molar 
concentrations  and  measured  surface  areas,  the  diffusion  coeffi¬ 
cient  for  the  LFP  cathode  and  LCO  anode  were  4.5  x  10-12  and 
8.3  x  10-11  cm2  s-1,  respectively.  These  diffusion  coefficient  values 
are  within  one  order  of  magnitude  as  those  reported  by  others 
[5,6,25,27,36,39,41-45]  reinforcing  the  validity  of  the  inference 
above. 


With  knowledge  of  the  diffusion  coefficient  it  may  be  possible  to 
estimate  the  average  particle  size  of  the  electrode  in  question  (as 
postulated  by  Dell'Era  and  Pasquali  [6  ).  Eq.  (14)  relates  the  effec¬ 
tive  particle  radius  (r)  to  the  Nernst  diffusion  length  (assuming 
spherical  particle),  <5,  provided  by  Srinivasan  and  Wang  [31  . 

where  the  Nernst  diffusion  length  is  related  to  the  diffusion  coef¬ 
ficient  by  the  definition  of  the  Warburg  resistance  (Eq.  (14)). 


RTS 

n2F2XSV2D 


(15) 


From  Eqs.  (14)  and  (15),  the  average  particle  radii  for  the  LFP 
cathode  and  LCO  anode  were  calculated  to  be  1.2  and  4.6  pm, 
respectively.  To  see  if  these  values  were  near  real  particle  sizes  for 
these  cells,  two  commercial  cells  (one  LFP  and  one  LCO  cell)  were 
disassembled,  their  anodes  and  cathodes  cleaned  in  EC  solvent, 
dried,  and  then  viewed  under  a  scanning  electron  microscope 
(SEM).  SEM  results  are  presented  in  Fig.  7. 

From  Fig.  7  it  can  be  seen  that  the  average  cathode  particle  radii  for 
the  LFP  cell  was  approximately  0.5-1.0  pm  while  the  average  particle 
radii  for  the  LCO  anode  was  approximately  12  pm.  In  the  case  of  the 
LixFeP04  cathode,  the  calculated  radius  is  almost  double  that  of  the 
SEM  result  while  the  calculated  carbon  anode  radius  was  about  one- 
third  of  the  SEM  result.  This  error  may  be  due  in  large  part  to  the 
orientation  of  the  electrode  materials  as  it  is  well  documented  that  the 
orientation  and  shape  of  electrode  particles  [46,47]  can  dramatically 
affect  the  diffusion  coefficient  [34,37,48,49].  Also,  in  the  case  of  the 
LiFeP04  cathode,  the  model  and  method  proposed  do  not  account  for 
any  carbon  painting  on  the  LiFeP04  particles,  commonly  performed  to 
increase  electronic  conductivity  [50  .  As  such,  while  it  may  still  be 
possible  to  calculate  the  particle  size  accurately  in  some  limited  cases; 
in  general  the  results  discussed  here  do  not  support  that  claim. 

The  electric  double  layer  capacitance  (ELDC)  observed  is  similar 
to  those  discussed  [51-53]  in  which  cathode  and  anode  interfacial 
capacitances  increase  as  their  potential  increases.  In  the  case  of  a 
full  cell,  as  the  cell  is  charged,  ions  leave  the  cathode  and  intercalate 
into  the  graphitic  anode  matrix,  corresponding  to  a  potential 
decrease  (Vanode  v.  Li-ion)  at  the  anode  and  a  potential  increase 
(Vcathode  vs.  Li-ion)  at  the  cathode.  From  the  data  provided  in  Table  4 
the  capacitance  values  designated  to  the  EDLC  remain  fairly  con¬ 
stant  across  the  cells  SOCs. 


Table  3 

Warburg  (diffusion)  parameters. 


Cell 

SOC 

Rw  (H) 

Tw  (s) 

Aw  (graphical)  (Q  s  0  5) 

Av  (calculated)  (Qs  0  5 ) 

DLi  (cm2  s’1) 

Average  particle  radius  (pm) 

LCO 

0.1 

0.055346 

62.1 

0.00700 

0.00702 

3.68  x  10^10 

4.60 

0.2 

0.050749 

47.5 

0.00810 

0.00735 

1.57  x  10  10 

0.3 

0.069428 

85.0 

0.00745 

0.00753 

8.69  x  10"11 

0.4 

0.07909 

91.5 

0.00740 

0.00826 

4.69  x  10-11 

0.5 

0.08157 

72.9 

0.00930 

0.00937 

2.57  x  10  11 

0.6 

0.09168 

100.0 

0.00995 

0.00900 

2.06  x  10  11 

0.7 

0.09900 

103.9 

0.00885 

0.00971 

1.36  x  10-11 

0.8 

0.12317 

182.4 

0.00770 

0.00912 

1.23  x  10  11 

0.9 

0.11849 

179.4 

0.00880 

0.00884 

1.06  x  10  11 

LFP 

0.1 

0.42096 

106.6 

0.03205 

0.04077 

1.07  x  1012 

1.20 

0.2 

0.23455 

143.1 

0.02010 

0.01960 

2.95  x  10~12 

0.3 

0.21352 

135.9 

0.01840 

0.01831 

3.83  x  10~12 

0.4 

0.24803 

127.8 

0.02035 

0.02194 

3.41  x  1012 

0.5 

0.24818 

145.7 

0.02045 

0.02056 

3.7  x  10  12 

0.6 

0.24692 

144.0 

0.01945 

0.02057 

4.49  x  10  12 

0.7 

0.24702 

110.0 

0.02145 

0.02355 

4.08  x  10  12 

0.8 

0.15503 

122.7 

0.01615 

0.01399 

7.98  x  1012 

0.9 

0.15229 

175.0 

0.01790 

0.01151 

7.25  x  10~12 
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Fig.  7.  SEM  imaging  of  (a)  LixFeP04-cathode  at  x  10,000,  (b)  Lii  *C6-anode  at  xlOOO 
and  (c)  Cross  section  of  same  anode  to  see  underneath  presumed  SEI  layer  at  x950. 


Table  4 

Double  layer  capacitance  as  a  function  of  SOC  and  surface  area. 


Cell 

SOC 

Cdi  (F) 

Calculated  anode 
specific  area  (cm2  g-1) 

Anode 
capacitance 
(uF  cm  2) 

Avg.  capacitance 
(uF  cm’2) 

LCO 

0.1 

0.0206 

237.45 

28.99 

28.47 

0.2 

0.0205 

28.88 

0.3 

0.0226 

31.76 

0.4 

0.0208 

29.23 

0.5 

0.0204 

28.78 

0.6 

0.0195 

27.41 

0.7 

0.0176 

24.80 

0.8 

0.0200 

28.11 

0.9 

0.0201 

28.37 

LFP 

0.1 

0.0081 

235.04 

14.76 

23.46 

0.2 

0.0108 

19.58 

0.3 

0.0152 

27.64 

0.4 

0.0198 

35.92 

0.5 

0.0123 

22.32 

0.6 

0.0117 

21.24 

0.7 

0.0118 

21.32 

0.8 

0.0120 

21.76 

0.9 

0.0147 

26.66 

area  and  volume,  along  with  the  measured  weight  of  the  cathodes 
and  anodes,  the  approximated  surface  area  could  be  calculated. 
Table  4  reflects  these  surface  areas.  Since  the  surface  areas  between 
LFP  and  LCO  cathodes  and  anodes  vary  by  several  orders  of 
magnitude,  the  overall  capacitance  of  the  cell  would  be  dominated 
by  the  smaller  anodes.  As  a  result,  the  capacitance  values  per  unit 
area  varied  as  a  function  of  SOC  from  0.292  to  0.274  F  nrT2  for  the 
LCO  anode,  and  0.147-0.266  F  m~2.  The  approximate  value  for  the 
EDLC  for  Lii_xC6  anodes  is  roughly  0.2  F  nrT2  has  been  reported 
[51-53].  By  this  comparison,  it  is  reasonable  to  conclude  that  it  is 
plausible  based  upon  a  rough  calculation  of  surface  area,  that  the 
anodes  of  both  cells  dominate  EDLC  as  described  by  EIS  and  PECM. 

5.  Conclusion 


By  using  a  physical  electric  circuit  model  (PECM)  that  accurately 
models  Li-ion  battery  performance,  characterization  of  the  electro¬ 
chemical  processes  occurring  inside  two  similar,  commercial,  Li-ion 
batteries  was  carried  out.  This  novel  method  of  applying  phenome¬ 
nological  or  physical-based  equivalent  circuit  modeling  to  EIS  data 
provided  unique  and  non-destructive  insight  into  the  electrochemical 
processes  inside  two  full-cell  Li-ion  batteries.  Both  cells  were  nearly 
identical  with  the  exception  of  their  cathodes  (LixFeP04  and  LixCo02). 
From  this  analysis  it  was  determined  that  the  cathodes  of  both  cells 
limited  the  exchange  current  densities,  the  LixFeP04  cathode  and 
Lii_xC6  anode  of  the  LFP  and  LCO  cells  each  limited  the  ionic  diffusion, 
and  the  anodes  of  both  cells  were  representative  of  the  double  layer 
capacitances  described.  It  was  shown  that  this  technique  provided 
this  information  which  may  be  useful  in  determining  the  limitation  of 
full-cell  batteries  without  disassembly. 
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To  validate  the  resulting  measured  EDLC,  the  surface  area  of  the 
electrodes  must  be  considered.  The  surface  area  was  estimated 
based  upon  the  average  particle  size  described  through  SEM  anal¬ 
ysis,  and  the  measured  weights  of  electrode  materials.  Through 
algebraic  manipulation  of  the  assumed  spherical  particle  surface 
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